Electron-phonon coupling and phonon self-energy in MgB2: 
do we really understand MgB2 Raman spectra ? 
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We consider a model Hamiltonian fitted on the ab-inttio band structure to describe tiie electron- 
phonon coupling between the electronic cr— bands and the phonon E2g mode in MgB2. The model 
allows for analytical calculations and numerical treatments using very large k-point grids. We 
calculate the phonon self-energy of the E25 mode along two high symmetry directions in the Brillouin 
zone. We demonstrate that the contribution of the cr bands to the Raman linewidth of the E2g mode 
via the electron-phonon coupling is zero. As a consequence the large resonance seen in Raman 
experiments cannot be interpreted as originated from the -E29 mode at P. We examine in details the 
effects of Fermi surface singularities in the phonon spectrum and linewidth and we determine the 
magnitude of finite temperature effects in the the phonon self-energy. From our findings we suggest 
several possible effects which might be responsible for the MgB2 Raman spectra. 

PACS numbers: 63.20.Kr, 63.20.Dj , 78.30.Er, 74.70.Ad 



I. INTRODUCTION 

The knowledge of the MgB2i electronic structure al- 
lows us to obtain a qualitative understanding of several 
peculiar features of this material . A crucial role is played 
by the a— band a-*"'" , formed by the in-plane boron-boron 
sp bonding. Due to the small interlayer coupling be- 
tween the boron layers, these bands have a two dimen- 
sional character and are weakly dispersing along the FA 
direction. Their corresponding Fermi surface sheets^ are 
two slightly warped cylinders, with axis perpendicular 
to the Boron layers. This peculiar topology results in a 
large contribution to the real and imaginary parts of the 
phonon self energy of the E2g phonon mode, an in plane 
displacement of the boron atoms. 

The large contribution to the real part of the phonon 
self-energy has spectacular consequences on the phonon 
spectrum: the phonon frequencies of the E2g modes un- 
dergo a reduction of roughly 20 meV along the FA direc- 
tions as predicted by ab-initio calculation a^i^i^ and mea- 
sured by high energy inelastic X-ray scattering-'-. Den- 
sity functional theory calculations of phonon^'^'^ spec- 
tra indicate that the softening of the E2g phonon fre- 
quencies when approaching the FA direction, even if 
strong in magnitude, is not as abrupt as would be 
expected^ for a pure two dimensional system having a 
Kohn anomaly. The experimental phonon dispersion^ 
along AL and FML^ confirms that the E2g phonon fre- 
quencies decrease gradually as the FA direction is ap- 
proached and the softening at momenta corresponding 
to the cylinders 2kp is very small. The Kohn anomaly 
might indeed be mitigated by the presence of a k^ band 
dispersion and a finite temperature. In this paper we 
investigate the magnitude of these two effects and dis- 
cuss their relevance in the interpretation of experimental 
data. 

A large contribution due to the electron-phonon cou- 
pling is also associated with the imaginary part of the E2g 
phonon self-energy, the phonon linewidth. Inelastic X- 



ray scattering experiments and theoretical calculations^^ 
show an anomalously large broadening (^ 20-30 meV 
FWMH) of the E2g mode along the FA direction only. 
According to ref>2, the broadening of this mode is almost 
temperature independent, but the spectra displayed in 
fig. 4 of ref.^ do not allow for a definitive conclusion 
since the J^2g rnode has a very small structure factor and 
it is seen only as a shoulder of the close Ei^ mode. 

Raman data show a completely different behaviour. 
Raman experiments probe excitations at small momen- 
tum transfer, close to the F point of the material. The 
maximum momentum transfer is qexp = '^lught where 
c[iight = ^ and A is the wavelength of the incident 
light. Most of the experiments are performed with a 
514.5 nm (2.41 eV) argon lascr^*' which corresponds to 
q^^p = 1.3 * lO-^ag-i « 0.002FM, (oo = 0.52921 is 
the Bohr radius). This region is inaccessible to X-ray 
measurement and as a consequence a direct experimental 
comparison between the two techniques cannot be per- 
formed. Nevertheless it is instructive to compare Raman 
spectra with the X-ray data as close as possible to the F 
point. 

The Raman linewidth of the Ezg modei2Jii2*iiiiiSii& 
shows a very strong temperature dependence since it is 
~ 20 meV (FWMH) at 40 K and reaches almost 40 meV 
at room temperature, a factor of two larger than the one 
detected in inelastic X-ray data along the FA direction. 
Since, according to the calculation performed in^, the 
anharmonic broadening at room temperature is 1.2 meV, 
it cannot be responsible nor of the large linewidth neither 
of its strong temperature dependence. 

An unexplored cause of such a large temperature de- 
pendence of the linewidth might be the electron-phonon 
coupling. The electron-phonon coupling contribution to 
the phonon linewidth is indeed temperature dependent, 
(see eq. 0] this work). Nevertheless the dependence on 
temperature is usually assumed to be negligible, but no 
detailed studies have been performed on the subject. 

In this work we carefully analyze all the approxima- 



tions involved in the calculation of the phonon linewidth 
due to the electron-phonon coupling. We analyze the 
temperature dependence of the phonon linewidth and 
the effects of neglecting the phonon frequency in Allen 
formula. It would be highly desirable to estimate the 
magnitude of these approximations using ab-initio calcu- 
lation, but the task is almost prohibitive. In actual ab- 
initio calculations^'^ a finite number of k-points is used 
together with a ~ 0.025 Rydberg {^ 3000i^) smearing of 
the Fermi surfacei^,. Physical effects involving tempera- 
ture difference between 40 and 300 K are basically invis- 
ible to the calculation, since grids having at least 1000 
times larger number of k-points would be needed^^. In 
this case, the calculated electron-phonon coupling con- 
tributions and its temperature dependence in the indi- 
cated region would be masked by computational details. 
The convergence of ab-initio calculations with the num- 
ber of symmetry-irreducible k-points is particularly rel- 
evant for MgB2<iS., since only the weak warping of the 
two cylindrical a bands Fermi surface sheets prevents the 
linewidth from diverging. Moreover, if there were effects 
such as anomalies in phonon spectra generated by the 
2fc_F singularities^, they could be detected only using a 
very large number of k-points in the phonon frequencies 
calculations. As a consequence the use of a too small k- 
points mesh might affect the calculation of both the real 
and imaginary part of the phonon self-energy. 

For these reasons, in this work we study the behaviour 
of the phonon self-energy of the Fj2g mode due to the 
electron-phonon interaction between this mode and the 
CT— bands using a model Hamiltonian. The Hamiltonian 
is composed by the two cr— bands coupled to an harmonic 
dispersionless E2g phonon mode. The form considered 
for the cr— bands is fitted from ab-initio calculations^ in 
the region close to the T point. The phonon frequency is 
that of the E2g at F. The model is illustrated in detail 
in sec. |nl together with the form of the phonon self- 
energy in its real (phonon shift) and imaginary (phonon 
linewidth) parts. The simplified form of the model al- 
lows to calculate analytically the linewidth as q ^ F 
along any high symmetry direction. Moreover it al- 
lows numerical calculations using grids of Nk — 300'^, 
symmetry-irreducible k-points in any point of the Bril- 
louin zone, which are enough to see temperature effects 
in the phonon linewidth due to the electron-phonon cou- 
pling. 

In sec. mil we calculate the phonon linewidth exactly 
in the limit q ^ F both in its intraband and in its cr — cr 
interband contributions. We consider two cases, (i) q 
along the FA direction (sec. IIII Al) and (ii) q along FM, 
or generally along any direction in the {k^, ky) plane (sec. 
IIIIBjI . since the bands and the considered coupling are 
isotropic in the {kx,ky) plane. We discuss the relevance 
of the results for the interpretation of Raman spectra. 

In sec. IIVI we follow rei^ and derive Allen formula 
starting from the phonon self-energy, paying particular 
attention at the approximations involved. Subsequently 
in sec. 0we numerically evaluate the phonon self-energy 
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FIG. 1: cr— bands Fermi surface cylinders with projection over 
the kz — tt/c and kz — planes. 



along the FA and FM directions for the case of a k- 
independent electron-phonon coupling. We estimate the 
effects of temperature and of cr — cr interband transitions. 
We evaluate the magnitude of the different approxima- 
tions involved in the derivation of Allen formula. The 
numerical results are also used as benchmarks to judge 
the reliability of preceding ab-initio calculations^ in what 
concerns the number of /c— points used in the simulations 
and the value of the smearing parameter. 

Finally we question the attribution of the 77 meV peak 
to the E2g mode and we suggest other interpretation of 
the Raman experiments. 



II. MODEL 

Following reference^, the structure of the a bands close 
to the Fermi energy can be expressed as: 



fkn — eo ^ 2t 



_ cos[kzC) -TZy 



(1) 



where the index n label the heavy/light hole bands. The 
holes masses are mi — 0.59 (heavy holes), m2 — 0.28 
(light holes). The average energy along FA, t± — 0.094 
eV, gives the dispersion of the bands. The top of the 
a bands is eo = 0.58 eV. Note that k are expressed in 
atomic units and —TZv in eV with TZy = 13.605. The 
bands are measured respect to the Fermi energy. The 
Fermi surface sheets identified by the bands in eq. Q] 
are two warped cylinders (see fig. ^. The radii in the 
k^ =0 plane are fcg^^'^ = 0.13 a^^ and fcjf^''^ = 0.09 a^\ 
The radn in the k^ = ±7r/c planes are k^'^ = 0.18 a,^^ 
andfct7p'^ = 0.126a~^ 

The contribution to the v phonon mode phonon self- 
energy due to the electron phonon coupling can be writ- 
ten as, 



ni,(q,Wqi,) = — Yl l.9kn,k+qml 



/] 



k+qrn 



Jkn 



Nk 



k,rn,n 



€k+q?7i ^kn 



^qi/ ■ 



IT] 



(2) 



where N^ is the number of k-points, the sum is over 
the Brillouin zone and /k„ are the Fermi distribu- 
tion functions. The matrix element is g^^^ k+qm ^ 
{'kn\dV/Suq^\'k + qm)/ y/luj^ where Uq^, is the amph- 
tude of the displacement of the phonon v of wavevector 
q, ujqi, its phonon frequency and V the electron-ion in- 
teracting potential--'^ . 

The real part of the phonon self-energy is 



rj 7VT / J iyk7i,k+qrn I 

k,r??..n 



/] 



k+qm 



/kn 



^k+qr? 



^krj 



, . (3) 

where P stands for the principal value. With this defi- 
nition, Aq express the renormalization of the harmonic 
phonon frequencies due to electron-phonon coupling ef- 
fects. 

The phonon linewidth (FWMH) is twice the imaginary 
part of ni,(q, Wqy) divided by Nk, as it can also be in- 
ferred from Fermi golden rule: 



Tqt' — pj / , l5kn,k-(-qml 
k,?n,n 

• (/kn — /k+qm) '5(ekH 



qm 



Ckri - t^qz.) (4) 



In the following two subsections we calculate the 
phonon linewidth in the limit q — + analytically. We 
use formula eq. 0] choosing q along two high symmetry 
directions in the Brillouin zone: (i) the out of plane FA 
directions and (ii) the in-plane FM direction. We show 
that in both cases the phonon linewidth vanishes in the 
q ^ limit. 



III. MGBa RAMAN LINEWIDTH. 

A. FA direction 

We choose q along the FA direction and consider the 
limit for q going to zero: 



f^k+qm 



II / TTl \ 

Ckri — - ~l]Tly + 2t±qcsm{k^c) (5) 



where k}. — k^ + k"^. Eq. © is correct at order 0{q^). 
Choosing Ljqi, = 65 meV (i.e. the harmonic E2g phonon 
frequency at F) and substituting in eq. ^ we get: 



47r 



Tqi' — j^ / , l5kn,k+qml (/k+qm /kn) 



km,n 



■5 



( — - -l\TZy - 2t^qcsm{kzc) + Uq^ 



(6) 



This sum can be divided in two different contributions, 
one coming from intraband transitions {m = n) and the 
other from interband transitions (m =/= n). 

If n = TO and in the limit of q ^ 0, in eq. 10 in 
order for the S function to be satisfied we must have that 



'^qi/ = 2t±qcsin{k2c). The momentum used in Raman 
experimen1ji£ q^xp = 1.3 * 10~^ao~^. Raman scattering 
then samples a sphere in momentum space centered at F 
and of radius qexp- It follows then that 



I2ii 



qe 



1.6 meV ^ LUc 



65 meV 



(7) 



using c = 6.653ao. Thus the (5-function condition in eq^l 
is never fulfilled in Raman experiments. The contribu- 
tion to the linewidth due to intraband transition is then 
exactly zero. The general fact that an optical phonon 
mode cannot couple with electrons at F as long as only 
intraband transitions are allowed has been already noted 
in the footnote number 18 of refi^S. 

Choosing a finite q along FA and using eq. Q we can 
determine the values of q in the Brillouin zone for which 
the intraband contribution is non zero, namely 



q>q 



.intra 



2i^ 



0.052 ao" 



O.ITA 



(8) 



perp^ 



In this estimate we have used the phonon frequency of the 
E2g at F, being the phonon branches fairly flat along FA 
and q -^ . We have also assumed the expansion © at 
order 0{q^) to be correct. We will show in sec. IV CI that 
this limit is indeed correct using numerical calculation. 

We then consider the interband contributions (to ^ 
n) and q ^ F. In order for the argument of the delta 
function in eq. ® to be satisfied we must have that 



2tj_qcsm{k;:c) — - - 1 7^y 



TO„ 



^qi/ 



(9) 



In the case n — 1 and m — 2, we have (^^ii — 1) > since 

' '^ 1712 ' 

nil > 7712. As a consequence the largest value of q for 
which the delta function is non zero is determined by the 
condition: 



\2t±qc\ < ujqjy 



(10) 



which leads to the same condition as in the intraband 
transition case, namely 



q > O.ITA 



(11) 



Then we consider the term with n — 2 and to = 1. 
In order for eq. ((HJ to give a finite linewidth at T = 
OK, the following two conditions must be simultaneously 
satisfied: (i) the states e^i are occupied and the states ek2 
empty (and vice versa) and (ii) the delta function in eq. |^ 
must be satisfied. Recalling the Fermi surface topology of 
the two a bands, the first condition means that the sum 
is limited to the region of space included between the 
two warped cylinders. This region is included between 
the two cylinders having axes along FA and radii A:™'" 
and k"^p^ respectively. On the other hand for the second 
condition to be fulfilled we must have 



> 



^ll 0.5257^y 

TO.2 2t±C 



2t±c 



(12) 



where we have substituted (^ ^ 1) = —0.525. The con- 
straint imposed by the two Fermi functions (condition 
(i)) allows to replace A:|| with A:™'" in the inequality: 



kl> 



(A:J,"'")2 0.5257^y lo^ 



m2 



2tic 



> 



2i I c 2t I c 



O.irA (13) 



From eq. 1131 we conclude that even the term with n — 2 
and m = 1 in eq. is zero for \q\ < O.IFA. A similar 
equation can be derived for the case eki > and ek+q2 < 
0, so that intraband transition give no contribution to 
the Raman linewidth. 

We have shown that both the a bands intraband and 
interband contributions to the phonon linewidth via the 
electron-phonon interaction are zero for q along FA and 



\q\ < O.IFA w 0.052 flg 



B. In-plane momenta. 



(14) 



We now chose q in the k^ , ky plane and consider the 
limit for q going to zero, we have that : 

|fc||+<7|2 |fc|||2 
Ek+qm = fkn 'R-U H 'R-V (15) 

where we have chosen q along the FM direction. The 
phonon linewidth becomes: 

Iqv = 'JT' / ^ IS'kn.k+qml (,/k+qm ~ Jkn) 



km,n 



|k||+q|' |k|||2 
TZy H TZy + uj^u 



(16) 



We neglect terms of order q^ . We first consider intraband 
transition only, [m — n). In this case, we obtain: 
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(17) 



where kx is along FM. In order for eq. ll7l to give a non- 
zero value for 7q^ the (S— function must be satisfied so 
that 2qkxTZy — rrimLu. The two Fermi functions limit 
the sum in the regions of space with (i) ekm < and 
ek+qm > and (ii) ek,n > and Ck-i-qm < 0. In case (i) 
the sum over k is limited to the region included by one 
of the two cylinders, depending on the value of the index 
m. The maximum k possible is fcmax = ^top'^ for m — 1 

(2) 

and fcmax = fci^t," for m = 2. We can then substitute 



.(1) 



.(2) 



fcmax and fcmax in the 5— function condition in eq. El to 
obtain qi = 7.8 x 10" 
92 = 5.3 X lO-^flg 1 a 



■^a^' « 0.0125FM for m = 1 and 
0.0085FM for m = 2. In case (ii) 
the sum is limited to the region of space outside one of the 
two cylinders and, since ek+qm < then q > fc|| — fc^p'^ for 



fc™ for TO = 2, with g > in both 



m — 1 and q > fc|| 

cases. Inserting g = fc|| — k^'^ or g = fcy 

(5— function condition and solving for kx one gets one gets 



fct'^'p'^ in the 



q[ = 7.4 X 10 



10 



-3„-l 



"o 



0.012FM for TO = 1 and q'2 = 5.1 x 



0.008FM for to = 2. Finally the intraband 



contribution vanishes completely for \q\ < \q\ < qi. 



which is factor of 



iRm{qi,q2,q[,q'2} = 5.1 x 10" 

4 larger than the exchanged momentum Qoxp. = 1-3 x 



10 



-3„-l 



in Raman scattering. 



Note that this conservative estimate have been ob- 
tained using the E2g phonon frequency at F. In the 
{kx,ky) plane the E2g phonon modes are not degener- 
ate and both have phonon frequencies which are larger 
than the value at F^. The use of a larger Wq^ would lead 
to the vanishing of the phonon linewidth at a larger value 
of q. 

Then we consider the interband case, (to, ^ n). We 
start considering n — 1 and m ~ 2. In order for the 
Fermi function difference in eq. 1161 to be non zero one of 
the following conditions must be satisfied: (i) e^i > and 
e/c+q2 < 0, (ii) Efei < and ek+q2 > 0. The first condition 



means that \q\ > fcg^^'^-fc™ « OMa^^ = 0.064FAf (the 



fc^ = plane is where the surfaces of the two cylinders are 
closer), i.e. no contribution to the linewidth for momenta 
smaller than 0.064FM. 

The second condition leads to fcii < k^'^^ and |kj| +q| > 



Thus we have: 



-TZy> 



(fco" 



TO2 



^ Uy > ^''^ 
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TO2 



TOl 



-TZy 



meaning that the S— function condition: 
TZy + ujq = Uy 

TO2 TOl 



(18) 



(19) 



is never satisfied. Then we consider the case n — 2, 
TO. = 1. The Fermi functions in eq. 1161 give the following 
two conditions: (i) ek2 < and ek+qi > 0, (ii) ek2 > 
and ek-i-qi < 0. In case (i) we have that fcy < fc™'" and 
|k 4- q| > fc™^'^ and we get the same result of n = 1, 
m = 2 case (i). In case (ii) we have fc|| > fc™*" and 
|k-|-q| < k}^^"". We have: 



fc2 

TO2 



ny> 



(fc™")2 
TO2 



jLmax 

TZy > ^^TZy 
mi 



(20) 



and similarly to the case with n = I and to, = 2, the 
condition 1 191 is never satisfied. 

In this subsection we have demonstrated that for q 
in the kx , ky plane the linewidth vanishes at small mo- 
menta, the intraband contribution vanishes for \q\ < 
0.008FM while the interband contribution vanishes for 
\q\ < 0.06FM. The phonon linewidth along FM vanishes 
for 



|g| <0.008FMw4ge 



(21) 



Conclusions 



(Aq^) for the mode v at point q one sees that 



In sections IIII Bl and IIII Al we have shown that the 
phonon hnewidth due to the electron-phonon couphng 
is zero in an eUipsoid centered at T and having axes 
(7|l = O.OOSaQ^ ~ 4:qexp in the {k^iky) plane and qj\_ = 



Ic^u 



0.052an 



40gexp along FA. Both axes are larger than 



the largest momenta accessible with Raman scattering, 
9exp = 1-3 X 10~^a^^. Along FA the minimal momen- 
tum giving a final linewidth is an order of magnitude 
larger than goxp.Thus if the Raman experiment is pre- 
pared with a geometry consistent with an exchanged mo- 
mentum along the FA direction one should indeed find a 
zero linewidth for the E2g mode. Although this geometry 
is currently employed in most of the Raman experiments 
in MgB2, it seems that a large linewidth (« AOmeV) 
is detected in the Raman data published up to now. 
We therefore conclude that the broad feature visible in 
these experiments cannot be associated to a pure E2g 
phonon excitations whose linewidth is determined by the 
electron-phonon coupling of the £'29 at Qoxp- In the final 
section of the paper we put forward possible explanations 
for the experimental spectra. 



IV. ALLEN FORMULA 

The linewidth 7qi, can be related to the electron- 
phonon coupling^*' via a simple approximations. Namely, 
at temperature such that kbT ^ uJc^u or in the case of a 
temperature independent ^c^v, using the J- function con- 
dition ^(ek+qm — ekn — Wqi/) in eq. Q) one can substitute 
in formula Q) 



/] 



k-t-qm 



/kn 



"qy 



W, 



V 



df 



(22) 



If the temperature dependence in equation Q is weak 
than the Fermi function can be considered a step func- 
tion, so that: 



l^y 



A-KLO 

1^ 



k.m.n 



1 5k, 



,k+qml '^(ekn)^(ek+qm - ^kn - t^) 



(23) 
If To is the highest temperature for which the substitu- 
tion of the derivative of the Fermi function with a Dirac 
(5— function in eq. |23|is still correct, then eq. \^\s valid in 
the range of temperatures such that fcfcTo > ki,T 3> oj^v 
If 7qj/ is temperature independent, then the condition is 
simply T < Tq. It is worth noting that in the limiting 
case of a very large phonon frequency it might occur that 
/cftTo < ^qf and formula 1231 might be never valid. Since 
in practice one has phonon frequencies which are of the 
order of 300K or more, the only real condition of ap- 
plicability of eq. 1231 is that jq^, has to be temperature 
independent. 

From the definitions^ of the electron-phonon coupling 



""- 2^iV(0X 



(24) 



which is Allen Formula^'^. Allen formula allows to ex- 
tract the electron phonon coupling from the measured 
linewidth under the assumption that anharmonic effects 
are negligible. For MgB2 this condition is fulfilled along 
the FA directionl. 

In actual calculations, it is customary to neglect the 
frequency dependence in the S function in ea. H23|l . ob- 
taining 



47rcL;, 



7, 



qi/ 



qi/ 



Nk 



k,'7n,n 



l5k«,k+qmP'5(ekrO'5(ek+qm) (25) 



This assumption is unjustified at q = and leads to 
the wrong behaviour at F. Thus formula H25|) cannot be 
used to explain finite temperature Raman experiments 
due (i) to its (wrong) behaviour at F and (ii) to the lack 
of temperature dependence. The correct behaviors are 
included in expression Q). 



V. NUMERICAL CALCULATIONS 



A. Model for the electron-phonon coupling matrix 
element. 



We consider a model composed of the two a bands in 
eq.(Q coupled to the phonons through a k-independent 
coupling. The electron phonon matrix element is: 
ffkm,k-i-qn = gSm,n + ag{^- Smn) , whcrc m, n run ovcr the 
two a bands and a determines the magnitude of the in- 
terband transitions (a = correspond to the case where 
interband transition are suppressed). We assume only 
one dispersionless phonon mode whose phonon frequency 
is determined from the calculated E2g phonon frequency 
at F^, namely Wq = 65meV = 754K. Along the FA di- 
rection, as confirmed by inelastic X-ray scattering datcii, 
this approximation is fairly correct for the E2g mode. 
Considering only one phonon mode, from now on we drop 
the index ly from the linewidth definitions. 



B. Technical details 

In the following subsections we calculate the real (eq. 
|3Jl and imaginary (eq. ^J parts of the phonon self-energy 
(eq. [21) in the k^ , ky plane and along the FA direction. 

In the calculations of the real part of the phonon self 
energy along FM we consider a finite temperature and we 
implement eq. |21with a rj smearing of 350/'ir. This smear- 
ing is necessary to calculate the principal value in eq. |21 
Thus we extract the real part at the end. This procedure 
gives a faster convergence as a function of Nk ■ The sums 
are performed using a grid of Nk — 300^ for the two di- 
mensional case with t± = and Nk — 300'^ for the three 



dimensional case with t± ^ 0. In both cases the grids are 
formed by N^. symmetry-irreducible A:— points, obtained 
from a mesh centered at F and randomly displaced from 
the origin. 

In the calculation of the imaginary part we replace the 
Dirac delta functions with Gaussians of width a, namely: 



6{x) 



(26) 



We compute then 7q, 7q and 7" (using equations 1^ 1231 
and 1251 respectively) as a function of ct on a given mesh 
of Nk symmetry-irreducible k-points. We then repeat the 
calculation on grids with always higher Nk (with N^ up 
to 300'^) in order to perform the limits of Nk -^ 00 and 
CT — > 0. In this way we obtain the continuum limit. The 
comparison between the results for the phonon linewidth 
obtained using the different formulas allows to judge the 
reliability of the different approximations in the calcula- 
tion of the phonon linewidth. 

In figure 121 (a) and (b) we show the convergence of the 
linewidths 7q (q = 0.6rA) as a function of the Gaussian 
smearing a (expressed in Kelvin). In (a) we used a Gaus- 
sian smearing, in (b) an hermitian-Gaussian smearing of 
order 1^^. As can be seen, the dependence on the smear- 
ing is fairly weak for the largest mesh {Nk = 27x 10^). In 
this case a value of a included between 50 and 100 K gives 
a result almost converged to the continuum limit. From 
now on we adopt this mesh and this range of a values to 
obtain the continuum limit of our 7 v.s. a curves. 

Note that previous ab-initio calculations of the phonon 
linewidth^ have been performed with mesh of 27000 
symmetry-irreducible k-points. For the case of an Her- 
mite Gaussian smearing, as can be seen in the picture, 
this would lead to an error in the estimation of the 
phonon linewidth of the order of ~ 5%. 
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FIG. 2: Phonon linewidth calculated at T = ZQQK using in 
formula 21 as a function of the Gaussian smearing [a), of the 
number Nk of k-points used in the sum. The arrow indicates 
the value that can be extracted from the largest mesh calcu- 
lation. In (a) a pure Gaussian smearing has been used, while 
in (b) an Hermite- Gaussian smearing of order 1 



C. Phonon self-energy due to the electron-phonon 
coupling 

1. Effect of the hand- dispersion along kz. 

In the {kx,ky) plane, the band structure of eq. ^ is 
composed of two bands each of them formed by a free 
electron like dispersion. As a consequence in the i_L — 0, 
a = and T — QK case (purely two dimensional with 
non interacting bands and at zero temperature) one ex- 
pects to find two singularities (one for each band) at 2kp 
in the imaginary part of the phonon self-energy and in 
the first derivative of the real part of the phonon self- 
energy. This is shown in fig. 13 (a) and (b) (dashed lines) 
for the real and imaginary part respectively. In the real 
part the singularities in the first derivative are seen as 
cusps at 2fci?. At T = the slope on the right of each 
cusp should be infinite. A finite slope is obtained as long 
as a finite non-zero temperature is usedSli (even for small 



temperatures the slope at 2kF is not vertical) . These sin- 
gularities are originated by the behaviour of the response 
function in two dimensions at T = and are smoothed 
out at finite temperature^. In three dimension {t± ^ 0) 
the singularities should disappea r^^^i^^ . The real part be- 
comes continuous with no cusps and in the imaginary 
part the singularities are replaced by smeared continu- 
ous peaks. The level of smearing is determined by the 
three dimensional character of the system, in our case by 
the strength of t± . Since the sigma bands in MgB2 have 
a small t±^ it is important to determine how far is the 
system from the two dimensional case. 

As can be seen in fig. |3 (a) and (b) (continuous line) 
the singularities are strongly affected even in the case of a 
small ij_ . Indeed the real part presents only very smeared 
cusps corresponding to the three-dimensional 2kp posi- 
tions. Similarly the imaginary part presents two smeared 
peaks at 2kp. Even the small tj_ considered in this pa- 
per is sufficient to basically eliminate the effects of the 
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FIG. 3: Real (a) and imaginary (b) part of the phonon self- 
energy of the E2g mode for f x = (dashed lines) and tx = 94 
meV (continuous line) calculated for q along the FM direction 
and at r = 407^'. Interband transition have been suppressed. 



which originates from the interband contribution. The 
interband contribution drops to zero at g « O.OGFM, as 
was predicted in section Fill Bl and as shown in the inset 
of fig. H 
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FIG. 4: Imaginary part of the phonon self-energy of the E2g 
mode calculated at T = 40/(' for q along the FM direction 
and at T = 40A' using in formula El and |1] for the real and 
imaginary part respectively, a = corresponds to the absence 
of interband transition. Inset: Intraband contribution to the 
E^g phonon linewidth. The dashed lines is g = 0.06FM, the 
limit derived analytically in sec. IIIIBI for the vanishing of 
intraband transition. 



Along the FA direction the interband transitions are 
completely negligible. This can be seen in fig. |7| where 
the two curves with a — \ and a = are indistinguish- 
able on the scale of the picture. 



two-dimensional ^kp singularities. 

We also study the behaviour of the phonon linewidth 
along the FA direction for fj^ ^ 0. Along this direc- 
tion, the phonon linewidth vanishes for |q| < O.IFA, as 
demonstrated in sec. IIII Al This is confirmed by the nu- 
merical calculations reported in fig. [7| (continuous line). 
The phonon linewidth increases monotonically approach- 
ing qo — O.IFA from larger momenta and becomes sin- 
gular for q -^ qo^, due to the behaviour caused by the 
(5— function in eq. 0] 



Z. Effect of the interband transitions between the a bands 

In this section we consider the effect of interband 
transitions, choosing a — 0,1 and t± = 9AmeV. The 
calculated imaginary part of the phonon self-energy at 
T = AOK is illustrated in fig. ^ Besides the 2fci? features 
found in the a = case we found several other features 



3. Temperature effects. 

Besides a finite dispersion along the kz axis, a second 
effect responsible for the smearing out of the singular 
features at 2kp is the finite temperature. In fig. [S] we 
show the phonon linewidth for momenta along FM for 
T = AOK and T = 300K. Overall there is a very weak 
dependence on temperature. Finite temperature effects 
(between 40 and 300 K) in the {kx,ky) plane are larger 
close to the 2kF singularities (see insect fig. 01 . Nev- 
ertheless, when compared to the value of the phonon- 
linewidth, temperature effects are fairly small and negli- 
gible in the calculations of the phonon linewidth. More- 
over, as can be seen clearly in the inset of fig. |S1 the 
singular behaviour of the two dimensional 2kp feature is 
completely lost. 

For the FA direction the effect is even smaller, indeed 
the results of the calculations at T = AOK and T = 300K 
are indistinguishable on the scale of fig. [7\ 
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FIG. 5: Imaginary part of the phonon self-energy of the E2g 
mode calculated at T = 300A' and T — 40K for a = 1 and q 
along the FM direction using eq. |1] 
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FIG. 6: Phonon linewidth of the E2g mode calculated at T = 
300K and a — 1 for q along the FM direction using eq. 0| ( 
7q),eq|ia(7q)andeq.|lll(7q). 



4- Allen formula 

In fig. Eland fig. [7|we compare the linewidth calculated 
using 7q (eq. ^, 7q (eq. ^^ and 7q (eq. EH) along the 
FA direction and FM directions respectively. 

In passing from the 7q to 7q we have assumed the 
linewidth to be temperature independent. In the preced- 
ing section (sec lVCSl) we have shown that this is indeed 
the case, so that we expect 7q ~ 7q almost everywhere in 
the Brillouin zone. This is what is seen in fig. Eland fig. 
[7\ These two pictures justify the use of Allen formula for 
MgBa. 

From 7q, 7^ is obtained by neglecting the phonon fre- 
quency in one of the two (5— functions of eq. 1231 This 
approximation leads to unpredictable results which, as a 
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FIG. 7: Phonon linewidth at T = 300K for q along the FA 
direction. The a — 1 case is overlapped to the a = 0, the 
contribution from intraband transition is very small. The 
linewidth vanishes for q < O.IFA (dotted line). 



consequence, must be investigated case by case, since the 
magnitude of the effects produced by this approximation 
crucially depends on details of the band structure close to 
the Fermi level and on the value of the phonon frequency. 
As shown in fig. jBl for q along the FM direction this 
approximation is fairly well justified. On the contrary 
along FA (see fig. [TJ 7q and 7° display two completely 
different behaviors. This is mainly due to the fact that 
7q is singular at F, while 7° is singular for q — *■ (see 
sec. nil A|) . Moreover, as we have shown in sec. 7q = 
for q < O.IFA. The proper behaviour is recovered in the 
region 0.3FA < \q\ < 0.5FA where we find that 7" « 7q. 



VI. CONCLUSIONS 

In this work we have studied the behaviour of the 
phonon self-energy of the E2g mode, both in its real and 
imaginary part. Our conclusions can be summarized in 
the following three points: 

1. Suppression of Fermi surface singularities in 
phonon dispersion and linewidth: two dimensional 
systems display 2k p singularities in the phonon 
spectrum and linewidth. Naively one would ex- 
pect MgB2 to be similar, being the band dispersion 
along the kz axis very small. On the contrary we 
have shown in sec. IV C II that even such a small t± 
strongly suppress the 2kp singularities, so that the 
phonon spectrum becomes rather smooth, and the 
singularities in the phonon linewidth are removed. 
An additional effect (see sec. IVC3|) is given by fi- 
nite temperature which at 300 K completely washes 
out any feature in the imaginary part of the phonon 



self- energy. 

2. Behaviour of the phonon linewidth for q — > F; 
We have shown that the phonon Hnewidth, both 
in its intraband and interband contributions van- 
ishes in an ellipsoid centered at F and having axes 
q\\ = 0.008FM in the {ka;,ky) plane and q± = O.IFA 
along the kz axis. The two values are larger than 
the Raman momentum qexp, namely q\\ ~ 4^qexp 
and q± ~ 40gcxp- This calculation demonstrates 
that the huge linewidth seen in Raman experiments 
cannot be attributed to the E2g mode. 

3. Temperature dependence in 7qi, and reliability 
of Allen formula ^q^ and of 7° ; the phonon 
linewidth is almost temperature independent in the 
T=0-300 K region. Small temperature effects are 
detected close to 2kp, but always less than some 
percent of the total linewidth. Since the phonon 
linewidth is basically temperature independent the 
use of Allen formula 7q^ is justified in the full 
Brillouin zone. On the contrary the approxima- 
tions customary employed in ab initio calculations 
of neglecting the phonon frequency in one of the 
(5— functions in 7q,y, obtaining 7q^ is not always jus- 
tified. Along the T — A directions the neglecting of 
the phonon frequency in the (5— function shifts the 
singularity at q « OAT — A to the F— point, leading 
to a completely wrong behaviour which affects all 
the region with q < 0.3FA. 



VII. CONSEQUENCES FOR THE 
INTERPRETATION OF RAMAN SPECTRA 

An immediate result of these three points is that the 
interpretation of the broad feature seen in Raman spectra 
at 77 nie\ii2*ii*i2iiiiiiiSii& as a phonon excitation due to 
the Ei2g mode at F is not correct. Indeed we have shown 
in this work that the huge temperature dependence of 
the Raman linewidth cannot be explained by a temper- 
ature effect in the electron-phonon contribution to the 
phonon linewidth. In a preceding work2^ we showed that 
the anharmonic contributions to the phonon linewidth 
has a weak temperature dependence. As a consequence 
the temperature dependence found in Raman data re- 
mains completely unexplained. Besides its temperature 
dependence, the value of the Raman linewidth at F is 
not consistent with the theoretical findings. Indeed we 
have demonstrated that electron-phonon contribution to 
the phonon linewidth is zero in an ellipsoid centered in F 
and larger than the Raman exchanged momentum. This 
is not at all the case for what concerns Raman spectra. 

There are additional considerations, concerning the po- 
sition of the 77 meV feature, which seem to indicate that 
it is very unlikely that it can be interpreted as due to a 
phonon excitation at F. The calculated harmonic phonon 
frequency of the F,2g mode at F is indeed 65 meV, a 
value 15% smaller than the experimental result. This 



has lead several groups to the conclusion that the dif- 
ference might be due to anharmonic effectsi^?^^^?. A 
careful determination^^ of the anharmonic phonon fre- 
quency shift, explicitly taking into account three and 
four phonon vertexes and the scattering between different 
phonon modes at different q-points in the Brillouin zone 
gives a fairly small value of this quantity at F (-1-5% of 
the harmonic phonon frequency, w 3.12 meV), clearly to 
small to justify the feature at 77 meV. This is confirmed 
by inelastic X-ray data of two independent groups^*^ 
showing phonon spectra in good agreement with the har- 
monic phonon frequencies, suggesting small anharmonic 
effects. Unfortunately inelastic X-ray scattering is not 
possible at the zone center, so that a direct comparison 
of the spectra cannot be performed. 

In what follows we analyze several hypothesis that can 
be made in order to reconcile, theory and inelastic X-ray 
data with experiments with Raman data. 

A possibility is that the 77 meV feature in Raman data 
could be ascribed to a single resonant process involving 
the E2g phonon mode coupled to electronic excitations. 
This would be consistent with the asymmetric shape of 
the peak, reminiscent of a Fano resonance^S. As a conse- 
quence the position of the peak would not correspond 
to the E2g mode phonon frequencies but it would be 
slightly shifted to lower frequencies. The temperature 
dependence of the linewidth might be different in this 
case. In this case it is interesting to study the peak po- 
sition as a function of the energy of the incident light. A 
study of the dependence of the spectrum from the wave- 
length of the incident light has been performed in refjii. 
It is shown that as the wavelength is changed, the peak 
energy position remains basically the same, even if the 
shape changes substantially. In the same work, from the 
study of the depolarization ratio between parallel and 
perpendicular orientations of the incident and emitted 
light, it is concluded that the symmetry of the excita- 
tion cannot be that of a single E2g mode, supporting the 
idea illustrated in this paper that the Raman does not 
measure the E2g phonon excitation at F. 

An alternative scenario is that the Raman peak might 
be due to excitation of phonons which are not at the F 
point. Such an excitation can be activated by (i) the pres- 
ence of defects such as Mg vacancies, (ii) multi-phonon 
scattering. A defect breaks translational symmetry and 
it makes possible to observe in Raman spectra phonon 
excitations at non-zero momenta. Indeed in a similar 
system such as defected Graphite, due to the almost two 
dimensional character of the electronic structure and a 
strong electron-phonon coupling^S*^, the phonon at the 
K-zone boundary has a very strong signal in Raman spec- 
tra (known as the D peak). However a defect activated 
peak cannot explain alone the strong temperature de- 
pendence of MgB2 Raman linewidth between 40 and 300 
K. This temperature range is not high enough to change 
the population of a phonon at 77 meV. The tempera- 
ture dependence might be explained by a multi-phonon 
process such as the absorption of an acoustic phonon and 
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emission of an optical phonon with opposite non-zero mo- 
menta. Multi-phonon scattering is also seen in graphite'^- 
and is responsible for the G' peak observed in Raman 
spectra. 



tions were performed at the IDRIS supercomputing cen- 
ter (project 031202). 
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